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Summary 

Three-dimensional mixing length models were 
developed for the analysis of flow field a small distance 
(less than 10 coolant hole diameters) downstream of 
coolant injection through a discrete hole at a 30“ angle 
into a crossflowing free stream. The experimental data, 
obtained with a hot wire sensor to probe the flow field, 
were used to determine the local mixing length. Mixing 
length models for the coolant exit region, separated flow 
region, and the reattachment flow region were inferred 
from the local mixing length distributions. To verify the 
mixing length modeling, conservation equations in a 
theoretical analysis were simplified and solved 
numerically, using the mixing length models to estimate 
the turbulence effects and to predict the velocity and 
temperature fields. With small injection rate and 
constant wall temperature, the numerical results of the 
crossflow streamwise velocity component and the surface 
heat-transfer rate were consistent with the hot-wire 
velocity measurements and with the results of the surface 
film cooling effectiveness study. With more detailed 
experiments, these mixing length models can be improved 
and used in theoretical analyses to study the effects of 
coolant temperature and turbulence on the discrete hole 
surface film cooling process. 


Introduction 

Coolant injection through discrete holes at an acute 
angle to the external flow is widely used for turbine blade 
surface cooling. The convective heat-transfer coefficient 
along the surface downstream of coolant injection is an 
important factor in determining film cooling effective- 
ness. Usually this heat-transfer coefficient is determined 
experimentally by measuring the local heat transfer rate; 
it is sometimes derived from the results of a two- 
dimensional boundary-layer flow analysis. Because of the 
complicated three-dimensional flow field associated with 
discrete-hole cooling, especially in the region near the 
coolant hole, two-dimensional boundary layer flow ana- 
lysis is inadequate to predict the local heat transfer rate. 
A three-dimensional computational method to predict 
this surface heat transfer rate is therefore needed. The 
technique of modeling the turbulent effects due to the 
mixing between the coolant flow and the external 
crossflow has been a major factor in the success of 
current computational methods. Thus, the objective of 
the present study was to develop three-dimensional 
mixing length models for use in the analysis of the flow 
field a small distance (less than 10 coolant hole diameters) 
downstream of a coolant hole. 

Recently, many theoretical analyses and experiments 
have been performed to study the turbulent mixin g of 


discrete-hole film cooling. Experiments (refs. 1 to 3) 
indicate that the coolant mixing layer turns toward the 
wall surface near the coolant exit and that the flow 
separates from and reattaches to the wall surface a small 
distance downstream of the coolant hole. Theories for 
predicting the flow process appear in references 4 to 6. 
The numerical technique applied in references 5 and 6 to 
predict the surface film cooling effectiveness was 
described in reference 7. An anisotropy coefficient was 
used in reference 6 to approximate the normal and lateral 
turbulent effects. Turbulent mixing length models, with 
linear mixing length distributions between pivot points, 
were used in reference 8 to compute numerically the flow 
properties immediately downstream of a two- 
dimensional tangential slot injection flow process. The 
present authors intend to develop similar types of mixing 
length model for the three-dimensional flow field 
considered here. 

In this paper three-dimensional mixing length models 
are presented for the coolant exit region, the separated 
flow region, and the reattachment flow region. These 
models were derived from the calculated local mixing 
length distributions obtained from the measured 
turbulence intensities. The resulting mixing length models 
are then used in a simplified theoretical analysis of the 
flow field to assess the effectiveness of these models. For 
this analysis, the governing conservation equations are 
solved numerically using an existing finite-difference 
scheme. The computed velocity and heat-transfer 
distributions are compared with the experimental results. 


Symbols 

Cp constant pressure specific heat 

D coolant tube diameter 

Eb anemometer dc voltage 

EbQ anemometer dc voltage at zero velocity 

Cb rms value of anemometer ac voltage 

Fi 2 empirical constant, =0.59 

K heat conductivity 

k turbulent kinetic energy 

/ local turbulent mixing length 

lb mixing length at u>d 

Iq mixing length at y = 0 

l\ mixing length at pivot point 1 

I 2 mixing length at pivot point 2 

P local generation rate of e 

Rp crossflow free-stream Reynolds number per meter 

St Stanton number, equation (15) 

Sto initial Stanton number 
T local mean temperature 


Tp free-stream mean temperature 

Tj coolant mean temperature 

Tyf wall surface temperature 

t fluctuating temperature 

U local mean velocity along x direction 
u fluctuating velocity in x direction 

Up crossflow free-stream mean velocity 
U j coolant mean velocity 

V local mean velocity along y direction 

V fluctuating velocity in y direction 

W local mean velocity in z direction 

w fluctuating velocity in z direction 

X longitudinal coordinate 

Xi longitudinal location where x=D 

X 2 longitudinal location where x=2D 

Xg longitudinal reattachment location 

y normal coordinate 

y\ y coordinate at pivot point 1 

y coordinate at pivot point 2 
y^ y coordinate at y = S 

z lateral coordinate 

Zc lateral coordinate at edge of coolant hole 
Zi Uj = constant for z£Zi 

Z 2 I unchanged for z^Z 2 

ai empirical constant, =0.94 

o !2 empirical constant, =0.26 

6 velocity boundary layer thickness 

e turbulent kinetic energy dissipation rate 

6 coolant injection angle 

X coolant injection mass flow rate, =pjU j/Pp^p 
fjL absolute viscosity 

p local density 

pp crossflow free-stream density 

p j coolant density 


Theoretical Analysis of the Flow Field 

The flow field of interest in the present study is a small 
region within 10 x/D of coolant injection. A schematic of 
this flow region is shown in figure 1 . The profiles of the 
mean velocity, mean temperature, and turbulent mixing 
length of the incoming crossflow are assumed to be two 
dimensional. The coolant is injected into the crossflow 
through a circular tube at a 30° angle to the wall. 

The mathematical model for an analysis of this type of 
discrete-hole cooling process is given in detail in reference 
5. The steady, three-dimensional turbulent flow field is 
described with continuity, momentum, and energy 


equations. That analysis is also used here. The set of the 
governing equations, the boundary conditions, and the 
closure assumptions, which solve analytically the velocity 
and the temperature fields, are given in the appendix. The 
x-direction gradients of the stresses and enthalpy flux and 
all the pressure gradient terms were neglected in the 
conservation equations of this study. This allows the 
effectiveness of the three-dimensional mixing length 
models to be assessed in a simple manner. 

The closure assumptions (eqs. (A7) to (A14)) require 
further information on the mixing length distribution 
before they can be used in the conservation equations. 
Thus, wind tunnel experiments with coolant air injection 
were performed to establish the turbulent mixing length 
models. 


Experiment and Data Analysis 

Details of the experimental apparatus were previously 
described in reference 9. The principal components 
included the tunnel, the test section plenum assembly, 
and a Hilsch tube for coolant air supply. The tunnel itself 
was made of clear plastic sections with a flow area section 
of 15 by 38 cm. The coolant injection hole (1 .27 cm diam) 
was located in the lower wall of the test section. The 
boundary layer thickness was 1.65 ^ ahead of the 
coolant hole with a crossflow velocity Upoi 22.9 m/sec. 
The mwn velocity U and the turbulence intensity, 
along the crossflow direction were measured 
with a single-sensor tungsten hot-wire probe. The axis of 
the sensor was parallel to the wall surface and 
perpendicular to the crossflow streamwise direction. To 


z/D 



y/D 
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take the measurements at locations across the boundary 
layer, the probe was traversed normally from the wall 
surface towards the external crossflow. 

Measurements 

The bridge dc voltages and the rms of the ac voltage 
were measured at x/D=0, 0.5, 1.0, 1.5, 2.0, 2.5, 3.5, 5.0, 
7.0, and 9.0 within the symmetric plane (z/D = 0). Similar 
measurements were taken at selected lateral locations, 
z/D = 0.25, 0.50, and 0.75, at x/D=3.5 and 5.0 (fig. 1). 
These data were collected at injection rates X of 0.22 and 
0.37. 

Mixing Length Computation 

The local turbulence intensity («2)i/Vc7 was computed 
using the following equation: 

= (1) 

which has derived from King’s law. This measured 
turbulence intensity was used to compute the local mixing 


length with the modified Reynolds stress transport 
equation developed by Hanjalic and Laimder (ref. 10): 

where and o!2 were empirical constants (see ^pendix). 
The normal and lateral gradients of the U velocity 
component were derived from the velocity profile 
measurements. 

Results of the Experimental Mixing Length 

Examples of the mixing length distributions, as 
computed from equations (1) and (2), are given in figures 
2 to 7. Figures 2 and 3 show the distributions over the 
coolant hole exit for the two injection rates, and figure 4 
shows the distributions downstream of the coolant hole. 
Lateral variations in the mixing length distributions are 
shown in figures 5 to 7 for different injection rates and 
longitudinal locations. For the present range of injection 
rates (X<0.4), these experimental mixing length 
distributions could be approximated with segmented 
linear functions of }> with pivot points ^i, ^2» and 



Turbulent mixing length, l/lg 

Figure 2. - Mixing length distributions over coolant exit region. 
Coolant injection mass flow rate, 0.22. 



Turbulent mixing length, t/lg 


Figure 3. - Variations of mixing length over coolant exit 
Coolant injection mass flow rate, 0.37. 
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Normal distance from wall surface, y/D 




Figure 4. - Mixing length distributions within plane of symmetry 
at z/D = 0. Coolant injection mass flow rate, 0. 22. 



0 .5 1.0 1.5 

Turbulent mixing length, 


Figure 5. - Lateral variations of mixing 
length at x/D = 3. 5. Coolant injection 
mass flow rate, 0. 22. 



0 .5 1.0 1.5 

Turbulent mixing length, Ift^ 


Figure 6. - Lateral variations of mixing 
length at x/D = 3. 5. Coolant injection 
mass flow rate, 0.37. 



0 .5 1.0 1.5 

Turbulent mixing length, 


Figure 7. - Lateral turbulent mixing length 
distributions at x/D = 5.0. Coolant Injec- 
tion mass flow rate. 0. 37. 
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However, the basic shape of the mixing length 
distribution changes as a function of longitudinal and 
lateral position. Coolant injection (XsO.4) did not 
disturb the mixing length distribution near the boundary- 
layer edge, and the edge of the boundary layer was taken 
as the outer pivot point. For the crossflow free- 
stream mixing Jength was 0.26 cm, with the following 
conditions: Up=22.9 m/sec, 7> = 300 K, and 
1.54x106. 

The mixing length distribution was different from that 
of the upstream undisturbed boundary layer only in a 
region close to the coolant exit where constant mixing 
length was found. (See figs. 2 to 4.) The local mixing 
length in the boundary layer increased, especially within 
the symmetric plane, after the coolant exit (fig. 4). The 
maximum turbulent mixing length li and the near-surface 
(y/D = 0.04) mixing length /j within the z/D = 0 plane are 
shown in Hgure 8 . A large variation in mixing length 
occurred in a region {x/D^5 for X = 0.37) immediately 
downstream of the coolant injection. After this region 
zero mixing length near the wall surface l\ and constant I 2 
( were found at the x/D locations within the limit of 
the present study. The lateral variation of the maximum 
mixing length I 2 is shown in figure 9. The experimental 
data are well represented by straight lines. 

Slope dl/dy, where 0sy^y2 at x/D = 3.5 and 5, was 
also determined from the experimental results and 
plotted (fig. 10). The effect of the coolant injection rate 
on dl/dy may be neglected for the cases considered. 


Three-Dimensional Mixing Length 
Modeling 

From the preceding experimental results, the present 
authors categorized the local turbulent mixing length 
distributions within the flow field of interest into three 
models corresponding to the coolant exit region, the 
separated flow region, and the reattachment flow region. 
The present models are derived from the comparisons of 
the normal direction mixing length distribution with the 
mixing length profile within the upstream boundary 
layer. This mixing length profile is approximated by three 
linear functions of y with pivot points, >^2 = 0-3 5, /2 = 0. 1 8 
cm and y = 6 , / 5 -O .26 cm (fig. 11(a)). Because of the 
insufficient data in the lateral direction, some assump- 
tions are made in the following modeling. 

Coolant Exit Region 

For the coolant exit region (fig. 1 1(a)) turbulent mixing 
length immediately over the coolant hole exit must be 
prescribed for this modeling (it was derived from the 
experimental data in the present study). The turbulent 
mixing length profile can then be determined as it is 


— O — Mixing length at pivot point 1 

— 0 — Mixing length at pivot point 2 



Longitudinal distance from coolant hole, x/D 


(a) Coolant injection mass flow rate. 0. 22. 

(b) Coolant injection mass flow rate, 0. 37. 

Figures. - Experimental results of I 2 and Y 2 within z = 0plane. 


Coolant injection 
mass flow rate, 

A 



0 .5 1.0 0 .5 1.0 


Lateral distance from symmetric plane, z/D 
Figure 9. - Lateral variations of mixing length at pivot point 2. 



Lateral distance from symmetric plane, z/D 


Figure 10. - Lateral variation of slope (dl/dy); 0 ^ y S y 2 - 
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Longitudinal distance, x Turbulent mixing length. I Lateral distance, z 


(a) Coolant exit region. 

(b) Separated flow region, 
(c) Reattachment flow region. 


Figure 11. - Turbulent mixing length models. 


presented in the figure. Mixing lengths /q and l\ were 
assumed to be constant along the z direction for 0 <z < zj 
and then to decrease linearly after z became greater than 
Z\. The mixing length Iq is set to zero at the wall surface 
for z '^Zc> h h the same as those of the upstream 
profile. 

Separated Flow Region 

The separated flow region (fig. 11(b)), is characterized 
by a large mixing length l\ near the surface and increasing 
li at different x/D. A linear variation in dl/dy for 

^y^y >2 as a function of x along the symmetric plane 
and a constant l\ for 0<y^yi are assumed. Mixing 
length /q = 0 is also imposed at ,y=0. 

Along the lateral direction l\ changes linearly from the 
centerline value to zero at Z 2 - Linear variation in bl/dy for 
yi <y <y 2 as a function of z is also assumed inO<z^Z 2 - 
The mixing length profile retains the undisturbed 
boundary layer profile for z ^Z 2 - 


Reattachment Flow Region 

For the reattachment region (fig. 11(c)) both /q and li 
are zero. Thus, within the plane of symmetry, the mixing 
length distribution is defined by I 2 and ^2 once they are 
given. A linear variation in bUby (between the pivot 
points y=0 and y =^ 2 ) as a function of z is assumed for 
0^z^Z2 to define the lateral mixing length distribution. 
The undisturbed boundary layer mixing length profile 
exists in z ^22- 

In the previous models assumptions have been made at 
Z=Z\ or Z2> and the properties l\, I 2 , and y 2 are also 
considered to be given at the pivot points. To completely 
define the mixing length models requires more detailed 
experimental results to establish empirical relations for 
these parameters. For the following computational 
verification of the present proposed models, linear 
approximations over the present experimental y 2 {x), 
l\{x), and l 2 {x) in the symmetric plane are used. The 
assumptions for the determination of z\ and Z 2 will be 
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described, with the boundary condition assumptions, in 
the following sections. 


Numerical Method and Assumptions 

The numerical procedures described in reference 11 
were used to write a set of finite difference equations 
corresponding to the governing conservation equations 
presented in the appendix. The finite difference equations 
were then solved numerically with a noniterative 
marching technique in the x direction and the standard 
tri-diagonal matrix algorithm to obtain the velocity and 
temperature distributions. 

The boundary conditions of the velocity and turbulent 
mixing length distributions at the coolant exit (y=^are 
shown in figure 12. Constant velocity components (JJ, V) 
within the ellipse (major ^is = 1.6D, minor axis = 0.8 D) 
enclosed by zi and linear U variation along the z direction 
for zi^z^Zc are assumed. Assumptions of these 
properties along the x-direction in the z = 0 plane are also 
given in this figure. An upstream disturbance distance of 
0.2 D along the z = 0 symmetric plane, with the undis- 
turbed two-dimensional mean velocity and temperature 
profiles in the y-z plane at this location (x= -0.2 D), is 
also assumed. These boundary conditions are not 
changed with different rates of coolant injection. 

An ellipse centered at x=xj, z = 0 is assumed for Z 2 for 
xsxj. Constant Z 2 from xi to x^ and the experimental 
results (fig. 5) are then used to determine Z 2 for x^Xg. 
Constant io the z^Z 2 is also assumed, and a step 
change in y 2 occurs at z=Z 2 - The size of the 
computational domain in the y-z plane was 1.5 U by 2 Z> 
and the grid point dimension was 35 by 61 for the present 


computations. However, a variable Ax step size (300 steps 
for X = 0.22 and 5(X) steps for X=0.37) with smaller Ax 
over the coolant exit region was used for x/Z)<10. To 
reduce the numerical diffusion and to avoid the use of 
fiiw grids near the wal^urface, a local maximum value of 
dU/dz = 500, and dUldy (fig. 13), through linear 
interpolation between the value in z = 0 plane and those 
from the power law proffle at zz, were imposed in the 
numerical process. The dUldy at the wall surface (y=0) 
was determined from the universal velocity profile 
assuming a sublayer existed there. 

To solve the kinetic energy dissipation rate equation, 
the turbulent kine^ energy at the coolant exit was 
assumed equal to uVoti (see appendix). A turbulence 
intensity of 0.07 (from th^xperimental measurements) 
was used to compute this The ratio of the turbulent 
kinetic energy to its rate of_dissipation at y = 0 was also 
estimated by k/e = l/{0.‘i dU/dy). 



Figure 13. - Assumptions of dU/dy in z = Opianefor 
numerical computation. 
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Zj Zg Z 2 z/D 
Lateral distance, z/D 

Figure 12. - Velocity and turbulent mixing length assumptions over coolant exit (y = 0). 
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Normal distance from wall surface, y/D 


Results and Discussion 

To verify the effectiveness of the present turbulent 
mixing length models, some numerical results are 
presented and compared with experiments. 

The Longitudinal Velocity Component 

The numerical results of this velocity component are 
plotted in figures 14 to 16 for \ = 0.22 and 0.37. The 
present measurements are also shown in these figures for 
comparison. 

For X = 0.22 (figs. 14 and 15) a small U component was 
computed near the coolant exit. Good agreement between 
numerical prediction and experimental results is found 
after the coolant exit. But the numerical approach failed 
to predict the near-surface velocity at ^/i? = 0.5, the edge 
of the coolant hole (fig. 15). For large injection rates 
(X = 0.37), good agreement between the numerical results 
and measurements within the symmetric plane is found 
only after the flow reattached to the surface (fig. 16). 


The discrepency between the numerical results and 
measurements with large injection rate may be due to the 
simplicity of the present mixing length modeling. The 
linear distribution between two pivot points is not 
satisfactory near downstream edge of the coolant hole 
with X=0.37. Constants ^2 and I 2 might not be adequate 
to describe the mixing length distribution over the 
coolant exit region. 

Analytically, the relation k/e = l/(0.3 dU/dy), which is 
valid within a local equilibrium turbxilence boundary 
layer, was used as the boundary condition at j'=0. The 
velocity gradient was then computed with the law of the 
wall, assuming a laminar sublayer existed over the w^ 
surface. For a large injection rate such as X=0.37, the U 
profile was significantly different from that in a fiat plate 
turbulent boundary layer. This k/e assumption might not 
be valid at high injection rate. 

The normal velocity gradient was also imposed with a 
mayimiim value to eliminate the grid dimension effect in 
this study. To release this constraint, a large number of 
grid points should be used. Smaller step size in the x 



Figure 14. - Comparisons of U components in the z-^0 plane. Coolant injection mass flow rale, 0. 22. 


x/D = 7.0 
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(a) Longitudinal distance, x/D = 3. 5. 

(b) Longitudinal distance, x/D ■= 5.0 


Figure 15. - Comparisons of U component in lateral plane. 
Coolant injection mass flow rate, 0.22. 


direction would also improve the numerical results, 
especially in the separated flow region. More extensive 
numerical computations might predict more accurately 
the properties in the near coolant hole flow field, 
especially for high coolant injection rate. 

Surface Heat Transfer Rate 

The numerical temperature profile was used to 
compute the local surface heat transfer rate. The second- 
order three-data-point-forward finite difference method 
was used to evaluate the temperature gradient at the wall 
surface. The local Stanton number then calculated 


s,-fr £<wa^)>-0 

PpUfCpi Tp- Tyy) 

From the initial profile Sto= 0.0016 was found. It is 
close to the value obtained from the flat plate boundary 




Figure 16. - Comparisons of U components in the z plane. Coolant injec- 
tion mass flow rate, 0. 37. 


layer equation (ref. 12) with the_ present crossflow 
conditions (7>= 300 K, = 278 K, C/p= 22.9 m/sec, and 
/lp= 1.54x106). 

The computational results of surface heat transfer 
rates in the z = 0 plane with X = 0.22 and 0.37 are shown 
in figure 17, 
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Figure 17. - Stanton number variations in z ■= 0 plane. 

Coolant to tree-stream temperature ratio, Tj/Tp ■ 0.963; 

wall surface to tree-stream temperature ratm, 0. 926. 

For \ = 0.22 significant reduction in the surface heat 
transfer rate downstream of the coolant hole was found. 
The maximum heat transfer rate was computed near the 
downstream edge of the coolant hole, and it reduced 
rapidly to a nearly constant value in the separated flow 
region. The surface heat transfer rate decreased again 
when the flow reattached to the wall surface and then 
increased to a constant value at locations far away from 
the coolant injection hole. 

For X = 0.37 a higher surface heat transfer rate 
(approximately the upstream undisturbed surface value) 
was computed immediately downstream of the coolant 
exit. This heat transfer rate reduced rapidly to a small 
value (nearly zero) after the flow reattached to the 
surface. A similar result has been measured from 
experiment (ref. 13). In many existing discrete-hole 
surface film cooling experiments, lower cooling 
effectiveness in the separated flow region with optimum 
cooling effectiveness in the reattachment flow region was 
reported at X»0.4. The difference between the present 
computational surface heat transfer rates along the plane 
of symmetry at X = 0.11 and 0.37 is consistent with those 
experimental film cooling effectiveness. 

The numerical results of the Stanton number in the 
lateral direction for X = 0.22 are shown in figure 18. The 
Stanton number increases rapidly at z/£)»l locations. 
The low surface heat transfer region (St < Sto) increased 
to twice the coolant hole diameter after the flow 
reattached to the surface. These phenomena are similar 
to the surface film cooling effectiveness described in 
references 3 and 9. A larger disturbed temperature field 
than the velocity field may occur at downstream 
locations. Using a small, constant y-z computational 
domain with the boundary conditions given at zlD= 1 (in 
the present computation) may induce the large lateral 
temperature gradient near z!D=\. Using a large 
^mputational domain may reduce this lateral 
temperature gradient near the edge of the domain and 
provide higher computational surface heat transfer rates. 



Figure 18. - Stanton number variations in lateral direction. 
Coolant injection mass flow rate, 0.22; coolant to free - 
stream temperature ratio, Tj/Tp, 0.963; wall surface to 
free-stream temperature ratio, T^^/Tp, 0.926. 


Concluding Remarks 

Three-dimensional mixing length models have been 
developed from the measurements of turbulence 
intensities within the flow field in a small region 
downstream of the discrete hole coolant injection at a 30* 
angle into a crossflow. Applying these models to a 
simplified three-dimensional theoretical and numerical 
analysis, the analysis predicted the longitudinal velocity 
distribution within the flow field with small coolant 
injection rates. With constant wall temperature, the 
numerical results showed a large reduction in the surface 
heat transfer rate downstream of the coolant hole. The 
distributions of the surface heat transfer rate were also 
consistant with the results of existing surface film cooling 
effectiveness studies. 

The results of the present study also indicate that, with 
the aid of the mixing length models, existing three- 
dimensional flow field theory may provide a method of 
computing the surface heat transfer rate resulting from 
discrete-hole film cooling, including the effects of 
turbulence and the temperature of the coolant. 

These mixing length models can be further improved 
with more elaborate experiments using a multiple, hot- 
wire sensor probe to measure the properties within the 
three-dimensional flow field. These additional data are 
necessary to completely define the mixing length models 
proposed here and to investigate the accuracy of the 
closure assumptions in the present analysis. 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, March 14, 1983 
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Appendix — Conservation Equations 


The Cartesian directions x, y, and z were alined with 
the crossflow streamwise, vertical, and lateral directions, 
respectively, and the correspqndii^ mean _velocity 
components were designated U, V, and W. The 
continuity, momentum, and enthalpy equations 
describing a steady, three-dimensional turbulent flow 
took the following forms, where, as discussed in the 
Theoretical Analysis of the Flow Field Section, the 
jr-direction gradients of stresses and enthalpy flux and all 
pressure gradient terms were neglected. 

^(p^+|;CoT0+^W = 0 (Al) 


9 / at/ _\ d / dU \ 





±(pUW)+±ifVW)+±(j,ir^ 


d r 

(dW 

dVT 

\ 9 , 

dW ~^\ 

Ty[^ 

A dy + 

dz \, 




(A4) 


and 




£ 

dy 


(£ djr_ 

\Cp dy 


/ K dT — 


) 


(A5) 


Boundary Conditions 

For the numerical solution of equations (Al) to (A5), 
the upstream undistributed two-dimensional mean 
velocity and mean temperature profiles must be specified. 


These profiles could be supplied by either assumption or 
measurement. For the present study the mean velocity 
profile was measured, and the mean temperature profile 
was computed using 


T-T^ _ U 

Tp-Ty, Jjp 


(A6) 


The following table lists other necessary boundary 


conditions: 

Boundary 

y/D 

u 

V 

W 

T 

Wall surface 

0 

0 

0 

0 

Ty. 

Coolant exit 

0 

Uj cos B 

Uj sin B 

0 

Tj 

Top surface (free 

00 

Up 

0 

0 

Tf 

stream) 

Symmetric place 

yID 

^K/^^ = 0 

BV/Bz = 0 

0 


(z=0) 

z, x^O 

y/D 

Initial profile 

0 

0 

Eq. (A6) 


Since the present study was directed at a small region 
near the coolant exit, the tube wall effect on the velocity 
distribution at the coolant exit (0^x/D^2 andy/D=0) 
was included. However, uniform coolant temperature 
was used. 


Closure Assumptions 

To solve the governing equations with the specified 
boundary conditions, the constitutive equations for the 
turbulence terms are required for the momentum and 
enthalpy equations. For the present study the following 
Reynolds stress transport equations developed for the 
thin shear flow analysis (ref. 14) are used in the 
momentum equations ((A2) to (A4)). 




“2 \^y / 
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-,.2fl2 + or,-2,2/3(/y 
«2 / 


(from ref. 5) are assumed in equation (A5). 

To find the rate of dissipation, the following transport 
equation was used: 


(A9) 

yyde T^e -^de 3 f/ 0.07 pA:2 \ 3e "J 


(AlO) 

dV/ 0.07 P^2\ del , Pe , _pe2 


(All) 

(A14) 

where cxi=u^/k, =(wi^^ + ww 2 )l/ 2 ^ 

^12 = ^12(“1»«2)- 

and 

with the local generation rate 


A set of the values, «i = 0.94, «2 = 0*26, and F \2 = 0.59, 
which are applicable to 90® corner boundary layer flow 
within a rectangular duct of arbitrary aspect ratio, has 
been established from experiment (reh 14). These 
constant values are used here to correlate the present 
crossflow streamline t/2 measurements. The y-direction 
distribution of the local mixing length, as it is computed 
from equation (A7), can be approximated with segments 
of linear y function. TTw experimental local mixing length 
is used to evaluate t/2, w2, uv, and uw terms in the 
momentum equations. 

The heat-transfer counterparts 


^ —dU —dU -jdV —,dW 

P= -UV -5 UW -r y2 w2-r— 

dy dz dy dz 


(A15) 


The e profile at the edges of the domain of interest was 
evaluated with the boundary layer assumptions, 

e = 0.3 k (dU/dY) 

Near the wall surface, 

e = 0.09075 k^.S/0Al/y 


vt=- 


3.2 e dy 


wt= — 


1 k—^dT 

-W2 — 

3.2 € dz 


(A12) 

(A13) 


(ref. 5), was used. The boundary conditions at y =0 (the 
wall surface) are e = 0 and k/e = constant and at y » oo are 
6 = ep and de/dy = 0. 

The corresponding boundary conditions at the coolant 
exit are given in the Numerical Method and Assumptions 
section. 
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